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We consider two issues in the DT model of quantum gravity. First, it is shown that the triangulation space 
for D > 3 is dominated by triangulations containing a single singular (D — 3)-simplex composed of vertices with 
divergent dual volumes. Second we study the ergodicity of current simulation algorithms. Results from runs 
conducted close to the phase transition of the four-dimensional theory are shown. We see no strong indications 
of ergodicity breaking in the simulation and our data support recent claims that the transition is most probably 
first order. Furthermore, we show that the critical properties of the system are determined by the dynamics of 
remnant singular vertices. 



1. Quantum Gravity 

The problem of reconciling the two theories 
of general relativity and quantum mechanics is 
perhaps one of the most important and difficult 
issues in theoretical physics. A variety of ap- 
proaches have been adopted see eg. but it is fair 
to say that only limited success has been achieved 
so far. In this article we shall concern ourselves 
with one such approach in which quantum me- 
chanical fluctuations of the (Euclidean) geometry 
are encoded in a path integral. 



Z{G) = J V[g] exp[-S(g,G}} 



(1) 



We can take the action to be the standard 
Einstein-Hilbert term S (g) = Yft~G I d A x^/g R 
with gravitational coupling G. A variety of diffi- 
culties are immediately observed; firstly, G carries 
dimensions of inverse mass squared so that the 
model is perturbatively nonrenormalizable. Per- 
haps worse S is unbounded from below so it is 
not even clear that a stable vacuum exists. Fur- 
thermore, even when the functional integral is re- 
stricted to some fixed topology it is not clear that 
a unique choice exists for the functional measure 
T> [g\. The approach we have been following is to 
replace this continuum integral with a finite lat- 
tice sum in which these problems can be evaded or 



at least regulated. This approach is termed DT- 
gravity. Some of the excitement about this model 
has stemmed from the initial observations [^) that 
the four-dimensional theory appeared to exhibit 
a phase transition at finite bare Go- This opened 
up the possibility of defining a continuum theory 
in the vicinity of this nonperturbative fixed point, 
perhaps in the spirit of Weinberg's work on 2 + e 
gravity @. 

2. Dynamical Triangulations 

The basic idea is that one can divide D- 
dimensional space into a collection of equal vol- 
ume elementary cells. The simplest choice for a 
cell is that it is a D-simplex — a set of (D + 1) 
fully interconnected points with edge length a 
taken to be some invariant cut-off. A piecewise 
linear approximation to the manifold is then con- 
structed by gluing these simplices together across 
their (D — l)-dimensional faces. Additional re- 
strictions are also commonly employed in order 
to ensure that the neighborhood of any vertex is 
homeomorphic to a D-ball. The local scalar cur- 
vature density is then proportional to the number 
of simplices common to a given {D — 2)-simplex 
which is bounded if the total volume is held fixed. 
Furthermore, its integral, the lattice action is es- 
sentially just the total number of nodes in the 
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triangulation. 

The partition function eqn [l] is then replaced 
by the lattice sum 



E 



exp 



KqNq — K4N4, 



(jy< - V) 

AV 



(2) 



In practice the bare cosmological constant K4 
is tuned to achieve a mean volume (JV4) = V 
and the remaining coupling ko plays the role of 
a bare Newton constant Go. The parameter AV 
controls the size of the volume fluctuations. 

There are many unknowns in this model. In 
this paper we consider just two — the structure 
of the measure over triangulations and the ergod- 
icity of the simulation algorithms currently in use. 
As we shall see there are some intimate connec- 
tions between these. 
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Figure 1. Singular local volumes for D=4,5 



3. Singular Structures 

Consider the triangulation space of the fixed 
volume V ensemble for D > 3. Let us define 
the local volume Ui of some (sub)i-simplex as the 
number of simplices which contain that i-simplex. 
Let us also call this simplex singular if u>i —> 00 
as V — > 00. 

Our results can be summarized in two conjec- 
tures: 

• Every triangulation contains exactly one 
singular (D — 3)-simplex composed of (D — 
2) singular vertices. 

• The local volume cj (D — 3) ~ aVi while 
the singular vertex local volumes luq ~ bV 

Fig [l] consists of two plots which illustrate these 
facts. The upper plot shows the volume of the 
singular (D — 3)-simplex for D = 4, 5 against the 
two-thirds power of the total volume V. Clearly 
the data lies on a straight line verifying the first 
half of the second conjecture. The lower plot 
shows that the singular vertex volumes vary lin- 
early with the total volume confirming the second 
part of this conjecture. Indeed, we have observed 
similar results in dimensions D = 6, 7 and can 
state that all the descendent subsimplices of the 



original singular (D — 3)-simplex have local vol- 
umes which diverge linearly with the total volume 
in contrast to the sublinear growth observed for 
the (D — 3)-simplex. Simple geometrical argu- 
ments are given in [Q for the two thirds power 
in the latter case. Indeed, the coefficient a can 
also be predicted in terms of b in this model and 
appears to fit the data rather well. 

The immediate question that arises is why 
there are so many triangulations possessing this 
singular structure. It is simple to understand 
this by introducing the concept of a local entropy. 
Consider the set of simplices which form the local 
volume u>i of some i-simplex and remove the com- 
mon (i + 1) vertices. The remaining vertices and 
their interconnections then form a triangulation 
of a dual (D — 1 — i) sphere. For example, a circle 
S 1 forms the boundary of the dual face to a link in 
three dimensions. The local entropy ascribed 
to that i-simplex is then defined to be equal to 
the number of ways of assembling the simplices 
in the local volume or equivalently the number of 
triangulations of this dual (D — 1 — i)-sphere. 

Si = Z D _ x _i (uji) (3) 

This function is known to increase at least ex- 
ponentially fast with volume for spheres with di- 
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mension greater than unity. Thus, z-simplices 
with (D - 1 - i) > 2 i.e. i < (D - 3) can in- 
crease their local entropy by growing their local 
volumes. It is then intuitively reasonable that a 
single singular (D — 3)-simplex ultimately domi- 
nates composed of a cascade of smaller dimension 
secondary simplices. 
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Figure 2. Correlation of singular vertex volume 
with total vertex number close to the transition 



This structure then forms a nonperturbative 
background in the crumpled phase of the model 
around which small fluctuations occur. Prelim- 
inary studies in four dimensions show that the 
volume of these singular vertices decreases mono- 
tonically as the coupling kq increases. This is 
eminently reasonable since forcing more vertices 
into the system by increasing the coupling en- 
hances the contribution of less singular triangu- 
lations (Singular ones have a minimal number of 
vertices for fixed volume). 

Close to the phase transition a gas of rem- 
nant singular vertices exists whose fluctuations 
are highly correlated with the fluctuations of the 
total number of vertices. This is illustrated in 
fig ^ which shows a perfect (anti)correlation be- 
tween the singular vertex volume loq and total ver- 



tex number N for V = 32K, D = 4, AV = 100 
and kq = 2.516 (the pseudocritical coupling). 
Thus the peak in the node susceptibility is driven 
by the fluctuations in the singular simplex struc- 
ture. The state with (No) ~ 6400 has a large 
mean extent and essentially no singular vertices 
- a local volume of approximately 1000 means 
that the vertex has merged with the background 
distribution of vertex local volumes and can be 
no longer distinguished as a singular vertex. Cor- 
respondingly the state with (No) ~ 6200 has a 
small extent and contains a variable number (at 
least one) of remnant singular vertices which are 
well separated from the background distribution. 
The critical system appears to tunnel repeatedly 
between these two states. We will return to this 
issue when we describe our experiments to check 
ergodicity. 

4. Ergodicity 

Current simulation algorithms typically target 
some volume V and allow for some (usually rather 
small) fluctuations of size AV around this tar- 
get value. It is known that for AV — * oo the 
Monte Carlo algorithm used in the simulation is 
ergodic. That is all triangulations may be reached 
by the application of the elementary moves. How- 
ever, the ergodic properties of the algorithm for 
finite AV are unknown. Indeed for AV^ < D 
the moves are known to be not ergodic. Thus it 
is important to check ergodicity carefully when 
studying the theory at finite AV. In the light 
of recent claims that the four-dimensional theory 
has only a first order transition we have studied 
that model close to criticality, using large lattices 
and for three different values of the fluctuation 
volume = 10, 100, 1000. A dependence of ex- 
pectation values on AV^ would signal a breaking 
of ergodicity. One way in which this could occur 
would be to suppose that the paths in the triangu- 
lation space between two volume V triangulations 
required intermediate triangulations with much 
larger volumes. Barriers of height B (V) might 
exist between states. Ergodicity would then be 
broken if 

AV < B(V) (4) 
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Figure 3. Time series for V = 32K, AV = 1000, 
k = 2.516 



Figure 4. Time series for V 
kq = 2.56 
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